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Abstract 

We have developed a simple method to construct local equilibrium function for one dimensional 
lattice Boltzmann method (LBM). This new method can make LBM model satisfy compressible 
flow with a flexible specific-heat ratio. Test cases, including the one dimensional Sod flow and one 
dimensional Lax flow are presented. Favorable results are obtained using proposed new method, 
indicating that the proposed method is potentially capable of constructing of the local equilibrium 
function for one dimensional LBM. 
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1. Introduction 

The lattice Boltzmann method is a mesoscopic-based approach for solving the fluid flow prob- 
lems at the macroscopic scales. In the LBM, the key issue is the determination of the relaxation 
time r and local equilibrium distribution function f( e i\ r represents the relaxation time that 
distribution function / relaxes to the local equilibrium function /( e ?). By Chapman- Enskog ex- 
pansion, relaxation time r is linked to the viscosity and the pressure that satisfy the Navier-Stokes 
equations. /( e< ?) is local equilibrium function that is constrained by macroscopic conditions. What 
kind of flow that can be simulated is determined by the form of the f( eqS > . In this paper we mainly 
discuss a simple method to construct local equilibrium function /( ei ?). 

In 1993, F. J. Alexander, S. Chen and J. D. Sterling ([1]) firstly to proposed a for 
compressible viscosity flow. In 1994, Chen, etc. (Q) developed Alexander's /( e? ) and obtained 
a model without nonlinearity deviations for recovering Navier-Stokes equations. In 2003, Watari 
and Tsuhara ([HI) pointed out that Chen's model can only be deployed for small viscosity and 
temperature range. They introduced a global coefficients to construct that make the model 
stable in a larger viscosity and temperature range. For the global coefficients in the Watari's 
model, more discrete velocity is required to make the /( ei? ) get correct form. 2007, Qu, etc. 
(@) proposed a circular local equilibrium function. In Qu's model, the circular local equilibrium 
function cannot satisfy high order macroscopic statistic conditions that make the model can only 
recover to Euler equations. 2008, Li, etc. (Q) developed Qu's circular equilibrium function. In 
Li's model, the /( e? ) satisfy all the macroscopic conditions which make correctly recovering of 
Navier-Stokes equations is possible. In the development history of /( e? ) in LBM, we can find that 
the macroscopic conditions are crucial for the constructing of the /( e? ) . 

In this work, based on the idea of assignment matrix to assign macroscopic conditions, we 
obtain a simple method to construct a local equilibrium distribution function for compressible 
flows. In our method, we do not consider the continuous form of local equilibrium function, and 
the discrete form of /( e? ) is directly derived by assigning the macroscopic conditions, so that f^ eq ^ 
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do not contain free parameters. For numerical simulation, the differential form of the lattice- 
Boltzmann equation is solved by the TVD scheme ((HI)- Flows with weak and strong shock waves 
were simulated successfully by the present model. 

The paper is organized as follows. In Section^ we introduce basic macroscopic conditions of 
LBM. In Section [3l the detail of our method to construct local equilibrium function is described. 
Firstly, we introduce assignment matrix to assign the non-energy conditions to construct f. ■ , 
then, we use fixed specific rest energy to make the f\ eq ^ satisfy energy conditions, Section[5]presents 
numerical results and discussions. In the last Section, we make some conclusions. 

2. Boltzmann-BGK equation for compressible flow and macroscopic conditions 

The standard dynamical theory of mesoscopic model is described by the Boltzmann equation 
as following 

g+£-v/-J. (i) 

where / is the distribution function ([3]). The £ is the particle velocity vector, and J represents 
the collision term, we can simplify the collision term J to the following formulation (0j): 

§£ + £-V/ = i(/^)-/), (2) 

where t is the relaxation time and f( eq "> is the local equilibrium distribution function. 
The moments of /( e,J ) are 

r f^d£ = p, (3a) 



f leq) Z a d£ = fm a , (3b) 

J f {eq] iaii}d£ = pu a up+pS a0 , (3c) 

f^ZaZpZ-ydg = pU a UpU 7 + p{u 7 S a fj + Uf}8 al + U a 6f3 7 ), (3d) 

j f [eq) (\\t\ 2 + Q)dt = P E, (3e) 

/ f {eq) (\m 2 + ot; a dti = (pE+p)u a , (3f) 

/ f {eq) {\\£,\ 2 +OU&d£ = (pE + 2p)u a u f3 +p(E + R g T)6 a0 , (3g) 

where the £ is the specific rest energy 0) and R g is gas constant. Consider the connection 

between mesoscopic variable and macroscopic variable: 

\p<? = P, (4) 

—pc 2 + p( = ps, (5) 

where D is the spatial dimension, c is peculiar speed and e is the internal energy. For the polytropic 
gas, we have 

C=[l-?(7-l)]e, (6) 



c = v/2(7-lK (7) 
where 7 is the specific heat capacity ratio. 
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Applying Chapman- Enskog expansion (Q) and moments (l3"al) - ([3g| to equations we can 
recover the macroscopic Navier-Stokes equations without body force: 

| + v>.) = o, 

' >(l>a) ' V • (pu®u) = -Vp + pAu + (/i + //)V(V ■ u), 



dt 
d( P E) 



V ■ [{pE+p)u\ = V • (AVT) + V ■ (it - IT), 



9t 

where II = /i[Vit + (Vit) T — jj(V • w)i] s /i is viscosity, // is the bulk viscosity and A is thermal 
conductivity. We have 

p = pr, A = (1 + ^)R g pr. (8) 



3. Construction of discretization local equilibrium function 

In above section, we known that the continuous /( e< ?) satisfy the macroscopic condition (j3a]) - 
|3g| ). For the lattice Boltzmann method, we need to known the discrete that deploy to the 

(9a) 



corresponding lattice microscopic velocity. The moments of discrete are 



(9b) 
(9c) 
(9d) 

(9e) 

(9f) 

(9g) 



By equations (p?a| - ( |9"g| ), the macroscopic conditions can be classified as non-energy statistic 
conditions and energy statistic conditions. The non-energy statistic conditions are 



^ iia^ili = pU a U fj +pS a/ 3i 

i 

E/i e9) (^i^i 2 + oe M = (^+pK, 



AT 



\ pUaUpUj +p(u 7 5 af 3 + up5 ai + u a Sp 7 ) J 



(10) 



From equation (|§a|) - ( |9"g| ), energy statistic conditions are translations of non-energy statistic 
conditions at some energy level. Therefore, f[ eq ^ should only conserve the non-energy statistic 
conditions. For the number of non-energy statistic conditions that we need to conserve, we choose 
corresponding number of microscopic £j, i.e., the number of lattice microscopic velocity £j is the 
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same as the non-energy statistic conditions. The square matrix that made up of moments of 
microscopic velocity is 



V 



1 



1 

£2 a 

6/3 



£la£l/3 ^2a£,2fl 



\ £iq£i^£i7 ^2a£,20£,2j 



\ 



1 



(ii) 



Define the V 1 as the assign matrix, where the V 1 is the inverse of V, we have discrete form 
ft q) ^ 



/ A (e9) \ 



(eg) 



jV- 



(12) 



9d|) . However, we know the f^ eq ^ satisfy non-energy conditions, we can not make the /f e< ^ satisfy 



By equation (fT2")l . we can simply find that the V * V 1 Sn satisfy equation (j9"aj) . (|9bl) . (l9"cl) and 
(I9dp . However, we know the f^ 9 ^ i 
the energy macroscopic conditions 

( 



St 



pE 
(pE + p)u 



\ 



V (pE + 2p)u a u p + p(E + R g T)6 aP J 



For thermal flow, £ is the specific rest energy. In equation (|5e|) . (|3fj) and ( |3"g| , £ is a microscopic 
variable and £ is unknown that make the total energy E can not be calculated from the moments of 
the distribution function Therefore, the introduction of a fixed Q to assign f^ eq ^ is necessary. 
The introduction of fixed Q can also solve the problem of j\ eq ^ which do not satisfy energy 
conditions. By Qu's work ([9]), linearly assigning / 4 onto two energy levels Ci and C2, Ci = 
and £2 > s max are the minimum and maximum stagnation energy in the whole flow field. We have 

,(eq) _ Aeq) & - C 
J iXi •* i t 1 

f( e l) _ f( e «)J>_ 

g}, ([5]) and (IT2"1) . we can easily verify that /j^y satisfy 



By equation (|9a| 



EE4?=* 
EE/S^ = ^> 

E E fiZ t<*& = PUaUp + pSafS, 

i i 

E E A.O = pUallpll-f + piUjdap + W^5 a7 + U a 5^) 



(13a) 

(13b) 
(13c) 
(13d) 
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EE4S^I^ + 0)=^ (13e) 

3 i 

EE4?4i*i 2+ g)& = (i3f) 
EE4?(^i 2 +o)^^ = (^+2 P k«/3 +p(£+fl fl T)^. (13g) 

3 i 

By this stage, entire f\ eq ^ of the lattice Boltzmann method for compressible flow with a flexible 
specific-heat ratio is established. 

4. Calculation scheme 

4-1- Time discrete 

For time discrete, we calculates advection term e*j • V/,(ff, cc)At as an explicit finite-difference 
form and the collision term —{fife, x) — fife, a;)) as a implicit finite-difference form, and we 
got 

/r +1 - f? + Ate- • v/r = M[ej? +i + (i - o)j?], (w) 

where x < t 4 < x + Ax, Jf = At(f^\ t ,x) - frfex)), Jf +1 = (f + At, x) - f t (t + At, x)). 

and represents the degree of implicity, in this paper we set 6 to be 0.5. 

Z.L Guo and T.S Zhao (@) introduced a new distribution function to remove the implicity of 
the equation (|14p . and the new distribution function is 

ff,-/ 4 + ^(/,:-/| e9) ), (15) 

where 7r = At/r. Applying this new distribution function to equation (|14j) . we can vanish the 
implicity of the collision term and got 

= -Atel ■ V/T + (1 - 7T + 7T0)tf + 7T(1 - 6)f\ eq) ' n , 



f n+l = _L_ {g n+l +7T0f (e qh n + ly 



(16) 



Space discrete 
For space discrete, we define 

f i (I,J) = f i (x I ,yj), f i (I-l,J) = f i (x I -Ax,yj), (17) 

where I and J are node indexes. The second-order TVD scheme is 

where Fi(I + h, J) is the numerical flux at the interface of xj + Ax/2, and is defined as: 

HI + \>J) = Ft {I + \>J) + Ft {I + \ J), (19) 

where 

Ft(I+\,J) = F+(J, J) + \mmmod (AF+(I + l,J),AF+(I-±J)), (20a) 

Ft{I+\,J)=F-(I + l,J)- iminmod (AFr (I + l,J),AFr(I+Z,J)), (20b) 

F+(I, J) = \(e la + \e ia \)fi(I, J), (20c) 

Ff(I,J) = \{e la - \e ia \)fi(I,J), (20d) 



and 

The function min mod( , ) in (gDafl and (f20b|) is the flux limiter ([l3|, [lit). 



A*f(J + i J) = Ft(I + 1, J) — ff (J, J). (20e) 
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5. Numerical examples 



In this section, we will applied new f- eqS> to Sod shock tube flow and Lax shock tube flow to test 
capability of this new method. In numerical computations, the dimensionless form is preferred. 
There are three independent reference variables for normalization, which are the reference density 
PO) reference length Lq, and reference internal energy eo- With the three reference variables, other 
reference variables and non-dimensional variables can be defined as 

Lq » t „ x p „ u e 

u Q = y/e ,to = — ,t=—,x = —,p=—,u = — , e = — . (21) 
u t L Q po u e 

For one dimensional LBM, we choose discrete velocity as e\ = 1, e2 = —1, e% = 2, = —2. By 
the assignment matrix, we have 



fx 



(eg) _ 



12 + ~6~ + ~6 3~~ ~ 3 

r2 „, „, 3 „ 2 



./„„■) / c u x c u x u x ^u x 2 

^ ) = -H-^ + T2-- + - + — "3 



- , '■- - , ] (^) 

(eq) f C C V'X V'X V'X 
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12 12 6 12 6 



r(eq) ( C Zi^; C 1l x 1i x 1i x 1 

/4 ' = 9 \ 8~ + I2~T2" + 1T + 12~6 

5.1. Sod shock tube 

First, we use the Sod shock tube to test the proposed model. Sod shock tube is a famous 
unsteady flow, which includes shock wave. The initial condition is given as 

(pL,u Ll e L ) = (1,0,2.5) (-0.5<f<0) 
(p R , u R , e R ) = (0.125,0,2) (0 < x < 0.5) 

In this case, we set f = 10~ 4 and the (2 — 4. The mesh size is taken as Ax = 1/201 and 
the time step size is chosen as At = f/4. Before the waves propagate to the two boundary ends, 
the distribution functions at the boundary can be set as the equilibrium distribution functions 
computed from the initial value of macroscopic variables. Figure [T] shows the computed density, 
velocity, pressure, and internal energy profiles (symbols) at i = 0.22, exact solutions (solid lines) 
are also displayed in this figure. Clearly, the present results agree excellently well with the exact 
solution. 

5.2. Lax shock tube 

The second test case is the Lax shock tube. The initial condition of the problem is given as 

(PL, u L , e L ) = (0.445,0.698,19.82) (-0.5 < x < 0) 
(p R ,UR,e R ) = (0.5,0,2.855) (0 < x < 0.5) 

We set f = 10~ 4 and C2 = 30. The mesh size and time step size are taken to be the same as those 
of the Sod shock tube problem. The computed density, velocity, pressure, and internal energy 
profiles (symbols) at t = 0.14 are shown and compared with the exact solution (solid lines) in 
Figure [21 Obviously, the present results are very accurate. 
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Figure 1: Velocity (left up), density (right up), pressure (left bottom) and internal energy (right 
bottom) profiles of Sod shock tube. 



6. Conclusions 

In this paper, a simple method to construct local equilibrium function for one dimensional 
LBM was proposed. To recover the fluid equation, we classify macroscopic conditions as non- 
energy statistic conditions and energy statistic conditions. For non-energy statistic conditions, we 
introduce assignment matrix to assign non-energy statistic conditions to obtain the local equilib- 
rium function f^ 9 ^. For energy statistic conditions, we use fixed specific rest energy, to assign 
]\ eqS> to different energy level that make the satisfy energy statistic conditions. Numerical 
experiments showed that compressible inviscid flows with weak and strong shock waves can be 
well simulated by the present model. 
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